Spin polarization of light atoms in jellium: Detailed electronic structures 
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We revisit the problem of the spontaneous magnetization of an sp impurity atom in a simple metal 
host. The main features of interest are: (i) Formation of the spherical spin density /charge density 
wave around the impurity; (ii) Considerable decrease in the size of the pseudoatom in the spin- 
polarized state as compared with the paramagnetic one, and (iii) Relevance of the electron affinity 
of the isolated atom to this spin polarization, which is clarified by tracing the transformation of the 
pseudoatom into an isolated negative ion in the low-density limit of the enveloping electron gas. 
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Interests in spintronics arc on the rise from both sci- 
entific and technological points of viewi^ Since devices 
in spintronics involve active control and manipulation of 
spin degrees of freedom in solid-state systems, it is abso- 
lutely necessary to have a deeper understanding of funda- 
mental interactions between electron spins and its solid- 
state environments. In view of this situation, we are in- 
terested in a composite system of an atom immersed into 
the otherwise homogeneous electron gas (EG). 

In an isolated atom, the ground state obeys the Hund's 
multiplicity rule that requires the highest spin configu- 
ration compatible with the Pauli's exclusion principle. 
Physically this rule is interpreted as the consequence of 
an effectively larger nuclear charge in a higher spin config- 
uration due essentially to the exchange effect. 3 Similarly 
in a uniform EG, the same effect favors spin polarization, 
bringing about the spontaneous spin-symmetry breaking 
or the spin-density-wave state which was proven to be 
the ground state at arbitrary electron densities within 
the Hartree-Fock (exchange only) approximation^^ The 
correlation effect, however, acts in the opposite direction 5 
and this effect is so strong in an EG as to lead eventu- 
ally to the paramagnetic ground state for the majority of 
metals. 

This paper deals with the composite system of an atom 
immersed into EG. Investigation of atoms embedded in 
the EG in both their paramagnetio§iL2i2ii2iAi and spin- 
polarizedi2iiii^iiSii& states has a long history. However, 
to the best of our knowledge, some important features 
of the electronic structure of the spontaneously spin- 
polarized states of this system have not been addressed 
so far. More specifically, they include: (i) Formation of 
the spherical combined spin density /charge density wave, 
which slowly decays with the distance from the impurity; 

(ii) Significant shrinkage of spin-polarized pseudoatoms 
as compared with their spin-neutral counterparts, and 

(iii) Demonstration of the way how the spin-polarized 
states of the impurities turn into those of the negative 



ions of the corresponding isolated atoms as the density 
of the enveloping EG tends to zero. The purpose of this 
work is to elucidate the above points. 

We are concerned with an impurity of the atomic num- 
ber Z (a pseudoatom) embedded into the otherwise ho- 
mogeneous EG at zero temperature characterized by its 
electron-density parameter r s = (3/47rrio) 1 / 3 , where uq 
is the uniform density of the EG in the absence of the im- 
purity. In the spin-density functional theory (SDFT), 17 
the Kohn-Sham equation is written in atomic units as 



-(l/2)A + t£"(r)] W(r)=eiW(r), 



(1) 



where the spin index a takes either f or |, e, and tpi jCr 
are, respectively, the energy level and the wave function 
of a Kohn-Sham electron orbital, v e J^ iy) given by 
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Z/r+ / [n(r') - n ]/|r - r'| dr' (2) 



is the effective potential, where n(r) = n^(r) + nj,(r) is 
the local electron density, u^ c ([ri|, njj; r) defined as 



(3) 



is the spin-dependent exchange and correlation (xc) po- 
tential with E xc [n^, njj being the total xc energy of the 
system, and v xc (no) is the spin- independent xc potential 
at the uniform electron density no. The spin densities 
are self-consistently determined as 



(4) 



The energy of a pseudoatom is the difference between 
the energies of the EG with and without the impurity: 
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where 8[ a (k) stands for the derivative of the phase-shift 
of the angular momentum I of the wave-function for 
a state in the continuous spectrum in the potential in 
Eq. ©. In Eq. JSJ) the first term represents the con- 
tribution from the bound states, the second term comes 
from the change in the density of continuum states due 
to the interaction with the impurity, while all the rest 
are ordinary (S)DFT contributions to the total energy^ 
regrouped to insure the convergence of integrals. 
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FIG. 1: Total energy of the spin-polarized (solid curves) and 
spin-neutral (dashed curves) states of the B, C, N, and O 
pseudoatoms versus the EG density parameter r s . 

We have solved Eqs. self-consistently for the 

atoms in the first two rows of the periodic table immersed 
into the EG of various densities. For H, He, Li, Be, F, and 
Ne pseudoatoms, spin-neutral ground states have been 
found in the EG density range of 3 < r s < 14. For B, C, 
N, and O pseudoatoms, on the other hand, we found spin- 
polarized ground states at the density of the EG lower 
than a certain threshold values, while the ground state 
was spin-neutral at higher EG densities. These conclu- 
sions agree with those of earlier studies^ 

In Fig. n we plot the total energy of Eq. JSJ of the 
spin-polarized and spin-neutral lowest-energy states of 
the B, C, N, and O pseudoatoms within the local spin- 
density approximation (LSDA) to the SDFT using the 
parametrization of the correlation energy of Ref. yjj. In 
all the four cases, below a definite threshold value of 
the EG density, which is different for different impurity 
atoms, the spin-polarized ground state has persistently 
lower total energy compared with its unpolarized coun- 



terpart. 

Our method of breaking the spin-symmetry was to 
start with imposing the occupancy of the 2p bound 
state with 3 electrons with spin up and less than 3 
electrons with spin down. Then we let the system re- 
lax self-consistently to its ground-state. No unoccupied 
bound states would remain upon the achievement of self- 
consistency: The 2p bound states we had had partially 
filled would disappear in the self-consistent potential for 
spin-down electrons. For spin-up electrons, depending 
upon the sort of the impurity atom and the density of 
EG, this state would either remain and then be filled 
with 3 electrons, or it would disappear as well. The net 
spin polarization would remain finite in either case. 

The results of the calculated spin densities for the car- 
bon atom in the EG of r s — 6 are shown in Fig. |3 to- 
gether with the total electron density of the polarized as 
well as the unpolarized system. We note that at larger 
distances from the center, the amplitude of the Friedel os- 
cillations of the total density in the spin-polarized state is 
significantly smaller than that in the neutral state, result- 
ing in the effectively more compact pseudoatom. The lat- 
ter finding is consistent with results for isolated atomsiSi 
The inset in Fig. [21 shows the local polarization 



C(r) = [n T (r) - nj. (r)]/[n T (r) + nj.(r)]. 



(6) 



The oscillating and slowly decaying local spin polar- 
ization around the impurity together with Friedel oscil- 
lations of the charge-density represent a spherical com- 
bined charge-density/spin-density wave. We determine 
the total electronic spin of the pseudoatom as 



S = (1/2) /[n T (r)-n|(r)]dr. 



(7) 



In Fig. |21 (left panel), the total spin of Eq. J7J is plotted 
against the electron-density parameter r s . We conclude 
that there exists a finite net spin excess or spontaneous 
magnetization of the impurity in the EG at electron den- 
sities below the threshold values. The net electronic spin 
of Eq. Q depends on both the atomic number of the 
impurity atom and the EG density, which finds itself in 
contrast with the result for the net charge of the impurity: 
Due to the full screening of a charge in the EG, which is 
closely related to the Friedel sum rule, the pseudoatom 
charge is 



Z = - [n T (r) + n L (r) - n ] dr, 



(8) 
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FIG. 2: Deviation of the density of electrons with spin up, spin 
down, and the total electron density from no/2, no/2, and no 
(dashed, dotted, and solid curves), respectively, around the 
C atom in EG of r s = 6. The dashed-dotted curve repre- 
sents the unpolarized calculation. The inset shows the local 
polarization of Eq. © . All curves are multiplied by 47rr 2 . 
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FIG. 3: Left: Spin of an impurity versus the EG density 
parameter r s . Solid lines are the fittings of the data with 
Eq. Right: The number of electrons in the sphere of 

radius R for pseudoatoms in EG (solid lines); Z + (R/r s ) 
(dashed lines). 



which is uniquely determined by the sort of the impurity. 

While at intermediate densities of the EG the total 
spin of a pseudoatom is governed by complicated many- 
body interactions within the impurity atom-EG system, 
the trend in a pseudoatom's spin at low densities (large 
r s ) has a clear qualitative interpretation. Because of the 
positive electron affinity (EA) of the B, C, and O iso- 
lated atoms (0.010, 0.046, and 0.054 a.u., respectively 21 ) 
the limiting case of these atoms immersed into the EG 
at zero EG density are the negative ions (NI) of the cor- 
responding atoms. According to the Hund's rule, the 
populations of the 2p orbital are with 2 electrons with 
spin up ( 3 P), 3 electrons with spin up ( 4 S), and 3 elec- 



trons with spin up and 2 electrons with spin down ( 2 P) 
for B _ , C — , and 0~ ions, respectively, corresponding to 
the total spin of 1, 3/2, and 1/2, respectively, which is 
clearly satisfied in Fig. at large r s . On the other hand, 
the NI of the N atom is unstable although long living 
(EA=— 0.003 a.uiSi), and the slow growth of the spin of 
this pseudoatom between 1 and 3/2 at large r s can be 
understood as the competition between the NI 3 P and 
atomic 4 S states. 

In the right panel in Fig. |3| the integrated number of 
electrons in a sphere of radius R are plotted versus the 
radius of the sphere for the EG of r s = 14. The plateaus 
in the case of B, C, and O close to the number of elec- 
trons of 6, 7, and 9, respectively, prove unambiguously 
the NI character of the corresponding states, while for 
N this number is between 7 and 8, inferring a state in- 
termediate between an atom and NI. The growth in the 
number of electrons to the right from plateaus is due to 
the electron density approaching the constant value of no 
at large distances from the center. This figure also shows 
that for a low-density EG, electrons extra to an atom 
or NI, whichever supported in the zero-density limit, are 
pushed away from the center leaving a region of nearly 
zero electron density between the atom/ion and the re- 
gion of nearly uniform EG, where the number of electrons 
in the sphere of radius R is approximately Z + (R/r s ) 3 
(dashed curves). 

For the period 1 and the rest of the period 2 atoms the 
same arguments lead to the spin-neutrality of the corre- 
sponding pseudoatoms: In the case of H, Li and F, which 
also have positive EA, the acquisition of an extra elec- 
tron completes the outer shells, causing the correspond- 
ing pseudoatoms to be spin-neutral in the low-density 
limit. For a different reason but to the same effect, in 
the case of He, Be and Ne, the spin-neutrality holds be- 
cause of the non-existence of their NI (even unstable ones 
with sufficiently long life-timeSi) , while the correspond- 
ing atomic states have zero spin. 

The steep fall in the spin of a pseudoatom near the 
critical point seen in Fig. [3] is suggestive of a phase tran- 
sition of the second order with the power dependence of 
the spin on r s near its critical value r sc 



S w a(r s - r sc Y 



(9) 



In Tabled the best fit values of the parameters in Eq. © 
are listed. These value strongly suggest that the expo- 
nent [i is universal and equal to 0.5. 13 
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TABLE I: Best fit parameters in Eq. ©. 

At the intermediate EG densities between the thresh- 
old value and zero, the total spin of a pseudoatom ob- 
tained via Eq. (JJJ is not, generally speaking, a multiple 
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FIG. 4: Sensitivity of the energy of the C impurity to the 
choice of the xc potential. Solid (dashed) curves refer to the 
polarized (neutral) states. The xc potentials used: LSDA of 
Ref. EH (PZ) and of Ref. E3 (GL), and GGA of Ref. |g3 
(PBE). 



of 1/2, as seen in Fig. EJ This fundamental difference 
between an isolated atom and the present pseudoatom is 
brought about by the contribution of the infinite number 
of delocalized electrons in the latter case^ 

In order to make connection to the earlier works as 
well as to test the sensitivity of our results to the choice 
of the xc potential, we have repeated the calculations for 
the C pseudoatom within LSDA using the parametriza- 
tion of the xc energy of Ref. and also beyond LSDA 
within the generalized gradient approximation (GGA) in 
its PBE version. 22 In the former case our results for the 
unpolarized states reproduce those of Ref. 0- As shown 



in Fig. 01 regardless of the choice of the xc potential, we 
have been able to obtain the spin polarized ground state 
of an atom embedded in the EG. While the total energies 
of the both polarized and unpolarized states arc shifted 
depending on a specific approximation, their difference 
(i.e., the stabilization energy of the polarized state) does 
not show considerable sensitivity to the choice of the xc 
potential. 

In conclusion, we have performed the spin-density 
functional calculation of the spin states of the period 
1 and 2 atoms embedded in electron gas. For H, He, 
Li, Be, F, and Ne pseudoatoms, we have obtained the 
spin-neutral ground states in a wide density range of the 
electron gas. On the contrary, for B, C, N, and O pseu- 
doatoms, there occurs a transition into the spin-polarized 
state at a critical density of the electron gas which de- 
pends on the atomic number of the impurity. Both results 
are in accord with earlier studies. In the spin-polarized 
state, the pseudoatom is found to be of a smaller ef- 
fective size compared with its spin- neutral counterpart, 
which is a feature in common with isolated atoms. We 
also observe a combined spherical spin-density/charge- 
density wave which manifests itself as the Friedel-like os- 
cillations. In the limit of the low density of the electron 
gas the electronic structure of a pseudoatom is found to 
converge to that of the negative ion of the correspond- 
ing isolated atom. The electronic structure of the spin- 
polarized state is largely different from that of the spin- 
neural one, which will certainly have impact on such ap- 
plications as the stopping power of metals for ions and 
the residual resistivity of alloys. 

VUN and CSK acknowledge support by the Korea Re- 
search Foundation by Grant No. KRF-2003-015-C00214. 



1 S. A. Wolf, D. D. Awschalom, R. A. Buhrmanb, 
J. M. Daughton, S. von Molnar, M. L. Roukes, A. Y. 
Chtchenkanova, and D. M. Treger, Science 294, 1488 
(2002). 

2 I. Zutic, J. Fabian, and S. D. Sarma, Rev. Mod. Phys. 76, 
323 (2004). 

3 R. J. Boyd, Nature 310, 480 (1984). 

4 A. W. Overhauser, Phys. Rev. Lett. 4, 462 (1960). 

5 A. W. Overhauser, Phys. Rev. 128, 1437 (1962). 

6 C. O. Almbladh, U. von Barth, Z. D. Popovic, and M. J. 
Stott, Phys. Rev. B 14, 2250 (1976). 

7 K. N0rskov and N. D. Lang, Phys. Rev. B 21, 2131 (1980). 

8 M. J. Stott and E. Zaremba, Phys. Rev. B 22, 1564 (1980). 

9 M. J. Puska, R. M. Nieminen, and M. Manninen, Phys. 
Rev. B 24, 3037 (1981). 

M. J. Puska and R. M. Nieminen, Phys. Rev. B 43, 12221 
(1991). 

1 C. Menchini, M. I. Trioni, and G. P. Brivio, Phys. Rev. B 
67, 035408 (2003). 

2 R. M. Nieminen and M. Puska, J. Phys. F: Met. Phys. 10, 
L123 (1980). 

3 N. Stefanou and N. Papanikolaou, J. Phys.: Condens. Mat- 



ter 3, 3777 (1991). 

14 N. Papanikolaou, N. Stefanou, R. Zeller, and P. H. Ded- 
erichs, Phys. Rev. B 46, 10858 (1992). 

15 N. Papanikolaou, N. Stefanou, R. Zeller, and P. H. Ded- 
erichs, Phys. Rev. Lett. 71, 629 (1993). 

16 P. Mavropoulos, N. Stefanou, and N. Papanikolaou, Phys. 
Rev. B 58, 1096 (1998). 

17 O. Gunnarsson and B. I. Lundqvist, Phys. Rev. B 13, 4274 
(1976). 

18 W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 

19 J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 

20 R. J. Boyd, Nature 250, 566 (1974). 

21 H. Hotop and W. C. Lineberger, J. Phys. Chem. Ref. Data 
14, 731 (1985). 

22 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

23 Because of the delocalized electrons, we can follow the con- 
vergence of a pseudoatom to NI within LDA, while for an 
isolated NI LDA fails (i 9 . The screening at large distances 
eliminates the difficulty of LDA for NI, i.e., an electron feel- 
ing the -1 charge if the self-interaction is not subtracted. 



